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Abstract 

A new set of three-dimensional visualisations of a large-scale direct 
numerical simulations (DNS) of a turbulent boundary layer is presented. 
The Reynolds number ranges from Ree = 180 to 4300, based on the 
momentum-loss thickness 6 and the free-stream velocity Uaa- The focus 
of the present fluid dynamics video is on analysing the coherent vortical 
structures in the boundary layer: It is clearly shown that the initial phases 
are dominated by coherent so-called hairpin vortices which are character- 
istic remainders of the laminar-turbulent transition at lower Reynolds 
numbers. At higher Re (say Reg > 2000), these structures are no longer 
seen as being dominant; the coherence is clearly lost, both in the near-wall 
region as well as in the outer layer of the boundary layer. Note, however, 
that large-scale streaks in the streamwise velocity, which have their peak 
energy at about half the boundary-layer thickness, are unambiguously 
observed. 

In addition to visualisation with classical three-dimensional isosur- 
faces, the video is also rendered using stereoscopic views using red-cyan 
anaglyphs. 



Simulation data: www.mech.kth.se/~pschlatt/DATAI 

Introduction Turbulent boundary layers constitute one of the basic building 
blocks for understanding turbulence, particularly relevant for industrial appli- 
cations. Although the geometry in technical but also geophysical applications 
is complicated and usually features curved surfaces, the flow case of a canon- 
ical boundary layer developing on a flat surface has emerged as an important 
setup for studying wall turbulence, both via experimental and numerical studies. 



However, only recently have spatially developing turbulent boundary layers be- 
come accessible via large-eddy simulations (LES) and even fully- resolved direct 
numerical simulations (DNS). The difficulties of such setups are mainly related 
to the specification of proper inflow conditions, the triggering of turbulence and 
a careful control of the free-stream pressure gradient. In addition, the numerical 
cost of such spatial simulations is high due to the long, wide and high domains 
necessary for the full development of all relevant turbulent scales. 

Simulation setup In this contribution we consider a canonical turbulent 
boundary layer under zero pressure gradient. More specifically, visualisations of 
the data sets obtained via large-scale DNS [51 [B] are shown, pertaining to a spa- 
tially growing boundary layer. The inflow is a laminar Blasius boundary layer, 
in which laminar-turbulent transition is triggered by a random volume force 
shortly downstream of the inflow. This so-called trip force, similar to a tripping 
strip in an experiment ^6J , is located at a low Reynolds number to allow the flow 
to develop over a long distance. Thus, the simulation covers a long, wide and 
high domain starting at Reg = 180 extending up to the (numerically high) value 
of Ree = 4300, based on the momentum thickness 9 and free-stream velocity 
Uoo- Fully developed turbulent flow is obtained starting from Reg « 500. The 
chosen numerical resolution for the fully spectral numerical method |2] is high 
enough to resolve the relevant flow structures; in the wall-parallel directions 
Aa;+ — 9 and Az+ = 4 is achieved. For the DNS, the whole simulation domain 
requires a total of 7.5- 10^ grid points in physical space (8192x513x768 spectral 
modes), and was thus run on a large parallel computer with 4096 processors. 

Flow characteristics Turbulence statistics obtained in the boundary layer 
such as mean profiles, fluctuations, two-point correlations, etc., of the flow com- 
puted by both DNS and LES are in good agreement with other simulations and 
experimental studies; see the comparisons presented in Refs. |31[51in]. When it 
comes to spectral information recorded in the boundary layer, footprints of two 
very distinct spatial (and temporal) scales are detected, i.e., the well-known 
turbulent streaks close to the wall scaling in inner (viscous) units, and long 
and wide structures scaling in outer units. These large-scale structures persist 
throughout the boundary layer from the free-stream down to the wall and act 
as a modulation of the near-wall streaks. It is thus clear that the interplay of 
these two structures, scaling in a different way, is strongly dependent on the 
Reynolds number, i.e., the downstream distance. 

Visualisation of coherent vortical structures However, the characteris- 
tics of these large-scale structures, and their relation to actual coherent struc- 
tures present in the flow are not entirely clear, in particular for higher Reynolds 
numbers. Recent studies, summarised by Adrian |T] and also observed in a DNS 
of low-i?e boundary layers [5| suggest a dominance of hairpin-shaped vortices 
of various sizes throughout the boundary layer, in accordance with the early 
suggestion by Theodorsen (1952). Hairpins are well-known structures in tran- 
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sitional flows, e.g., appearing as the results of shear-layer roll- up. It is now 
interesting to see whether hairpin-like structures persist in fully turbulent flow 
(either in the outer region or close to the wall), or whether they are mainly 
restricted to low-i?e and/or transitional flows. Thus, we would like to study the 
structural characterstics of a turbulent boundary layer as a function of Reynolds 
number, starting from the transitional phase up to fully turbulent flow at higher 
Re than available in previous studies. 

The videos in this contribution demonstrate the vortical structures arising 
in the boundary layer as function of the Reynolds number, Reg, i.e., the down- 
stream distance. In all frames isocontours of negative A2 [3] identifying vortical 
structures are coloured by the wall distance. The video consists of two parts, 
essentially showing a trajectory of the camera through an instantaneous turbu- 
lent field. The first part shows classical isosurfaces at — —0.05, coloured by 
the wall distance. The second part, based on the same trajectory, is rendered 
using a stereoscopic view with red-cyan anaglyphic images. To view properly, 
standard red-cyan glasses are required. 

Physically, laminar-turbulent transition is induced by trip forcing close to 
the inlet in a similar fashion as in experiments. The subsequent breakdown to 
turbulent flow is characterised by the appearance of velocity streaks and un- 
ambiguous hairpin vortices, which are seen to dominate the whole span of the 
flow, see e.g., Ref. [T]. The hairpin vortices increase in number, and individual 
distinctive heads of such vortices are clearly visible for some distance down- 
stream at Reg ~ 800. This feature of low-i?e turbulent flow is also put forward 
in Ref. [8j, there denoted as "forest of hairpins". We can thus confirm that 
at least at low-i?e hairpins are indeed the dominant structure in a turbulent 
boundary layer. However, as the Reynolds number is further increased above 
about 1300, the scale separation between inner and outer units is getting larger, 
and the flow is less and less dominated by these transitional flow structures. At 
Reg = 2500 only isolated instances of arches belonging to hairpin vortices can 
still be observed riding on top of the emerging outer-layer streaky structures. 
But the dominance of hairpin-likc structures is clearly lower than in the previ- 
ous figures closer to transition. This effect is even increased by considering the 
highest present Reynolds number, Reg = 4300. Then, individual hairpins or 
arches cannot be seen any longer. The boundary layer is now truly turbulent, 
and the outer layer is dominated by large-scale streaky organisation of the tur- 
bulent vortices. Ongoing analysis of the flow structures close to the wall (i.e., 
the near-wall cycle) reveals charactersitic oscillations of the near-wall streak, 
mainly in a sinuous manner leading to a staggered appearance of the vortex 
cores [7]. 

Acknowledgements Computer time was provided by the Swedish National 
Infrastructure for Computing (SNIC) with a generous grant by the Knut and 
Alice Wallenberg Foundation (KAW). 



3 



References 

[1] R. J. Adrian. Hairpin vortex organization in wall turbulence. Phys. Fluids, 
19(041301):1-16, 2007. 

[2] M. Chevalier, P. Schlatter, A. Lundbladh, and D. S. Henningson. SIMSON - A 
Pseudo-Spectral Solver for Incompressible Boundary Layer Flows. Technical 
Report TRITA-MEK 2007:07, KTH Mechanics, Stockholm, Sweden, 2007. 

[3] J. Jeong and F. Hussain. On the identification of a vortex. J. Fluid Mech., 
285:69-94, 1995. 

[4] P. Schlatter, Q. Li, G. Brethouwer, A. V. Johansson, and D. S. Henningson. 
Simulations of spatially evolving turbulent bormdary layers up to Reo = 
4300. Int. J. Heat Fluid Flow, 31:251 261, 2010. 

[5] P. Schlatter and R. Orlii. Assessment of direct numerical simulation data of 
turbulent boundary layers. J. Fluid Mech., 659:116-126, 2010. 

[6] P. Schlatter, R. Orlii, Q. Li, G. Brethouwer, J. H. M. Fransson, A. V. Jo- 
hansson, P. H. Alfrcdsson, and D. S. Henningson. Turbulent boimdary layers 
up to Ree = 2500 studied through numerical simulation and experiments. 
Phys. Fluids, 21(051702):l-4, 2009. 

[7] W. Schoppa and F. Hussain. Coherent structure generation in near-wall 
turbulence. J. Fluid Mech., 453:57-108, 2002. 

[8] X. Wu and P. Moin. Direct numerical simulation of turbulence in a nominally 
zcro-pressure-gradient flat-plate boundary layer. J. Fluid Mech., 630:5-41, 
2009. 



4 



